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We present the fully differential production and decay of a W' boson, with arbitrary vector 
and axial- vector couplings, to any final state at next-to-leading order in QCD. We demonstrate 
a complete factorization of couplings at next-to-leading order in both the partial width of the W' 
boson, and in the full two-to-two cross section. We provide numerical predictions for the contribution 
of a W' boson to single-top-quark production, and separate results based on whether the mass of 
the right-handed neutrino ur is light enough for the leptonic decay channel to be open. The single- 
top-quark analysis will allow for an improved direct W' mass limit of 525-550 GeV using data from 
run I of the Fermilab Tevatron. We propose a modified tolerance method for estimating parton 
" distribution function uncertainties in cross sections. 
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^ , I. INTRODUCTION 

Since the introduction of the standard model, many extensions have been proposed that involve enhanced gauge 
' ' ' symmetries. One common feature of these models is the prediction of additional gauge bosons, generically called 
W and Z' bosons. In non-universal and top-flavor models, the W gauge boson arises from a new SU(2)l sector that 
distinguishes between generations of fermions |^-|^ , and may treat quarks and leptons differently |^-^ . The W boson 
could be the lowest Kaluza-Klein mode of the W boson [Q, or a heavy mass eigenstate in non-commuting extended 
<^ ' technicolor jsj. In a left-right symmetric model the W' and standard- model W bosons are remnants of a broken 
' SU(2)^ X SU(2)^ symmetry [pHl2[. The resulting left-right mixing may be naturally suppressed by orbifold breaking 
of the left-right symmetry [^3|| , or by supersymmetric interactions . While indirect bounds may be placed on the 
masses and couplings of these various W bosons within the context of their explicit theories, it is always advantageous 
to search for these new particles directly. In order to facilitate a direct experimental search, and to achieve the most 
general contact with theory, we calculate the fully differential next-to-leading order cross section for the production 
and decay of a W boson to any final state with arbitrary vector and axial- vector couplings. 

The most general Lorentz invariant Lagrangian describing the coupling of a W' to fermions may be written as 

^ - jf MChP^. + Ctf^PL)W'f, + H.c. , (1) 

" where Pu.l = (1 i 75)/2, g is the standard model SU(2)i coupling, and the C"^^ are arbitrary couplings that differ 

for quarks and leptons. For a standard model W boson, — 0, and is either the Cabibbo-Kobayashi-Maskawa 
• ^ matrix or diagonal, for quarks or leptons, respectively. 

Most experimental searches have concentrated on W' bosons that decay via a purely left- or right-handed current. 
M In order to make contact with these results, we rewrite Eq. in the notation typical of left-right symmetric models 
. -"^ . as H 

C - ^7,7m {9Re'^ cos C Pr + Ql sin ( F^t/, Pl)W' f, + H.c. , (2) 

where C is a left-right mixing angle, and is a CP-violating phase that can be absorbed into V^. In this notation, 
gR(L) are the right (left) gauge couplings, and V^j:^ are generalized Cabibbo-Kobayashi-Maskawa (GCKM) matrices. 
In models where the W and W mix, the mixing angle ( is usually constrained to be small (|C| < a few xlO~^-10~^ 
1^). Hence, the search for purely right- or left-handed states appears well motivated. 

Many direct searches for W' bosons have been performed at hadron colliders. Most experimental analyses have 
considered left- or right-handed W bosons, with standard-model- like couplings, that decay into leptonic final states. 
Left-handed bosons, or right-handed W'j^ bosons in which the decay into a right-handed neutrino i/ji is kinemati- 
cally allowed, are constrained to have masses mw' > 786 GeV ]l7|-|2^ . If m^,^ > m^, the decay to vr is not allowed, 
and the right-handed bosons are only directly constrained by peak searches in the dijet data. Unless the W has 
greatly enhanced couplings to light quarks, the dijet data are limited by QCD backgrounds to providing a mass limit 
of 420 GeV (2|-|2^] 
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The only final state not measured by one of the experiments listed above involves the decay of a W' to a single top 
quark. In Refs. [p8p9|] it was pointed out that a deviation of measured the single-top-quark cross section from the 
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standard model prediction could be evidence of a new gauge interaction. In the context of models with W bosons, 
we wish to take this one step further, and propose that the experiments search for an explicit W' mass peak in the 
s-channel single-top-quark sample. We examine just how large the cross section into single top quarks can be at 
the Fermilab Tevatron and Large Hadron Collider (LHC) as a function of W mass, and determine how enhanced or 
suppressed couplings enter into the measurable cross section after cuts. We also determine the next-to-leading order 
distributions of the final-state jets to estimate the effect on the reconstruction of the W' mass peak. 

A W boson that propagates in the s-channel appears very similar to DrcU-Yan production. We might be tempted 
to use the next-to-next-to-leading order (NNLO) predictions for Drell-Yan (or the updates in Ref. However, 
these calculations are not adequate to predict the cross section for W' bosons which decay to quarks. There are 
two reasons for this. The first problem is that the final state effects are very large. Not only are there large QCD 
corrections to the quark final state, but the top-quark mass has a large effect on the branching fraction when it is in 
the final state. The second reason is that there are initial- and final-state interference terms that arise at NNLO that 
are of an unknown size. (Pure QCD processes also interfere at NNLO, e.g. with ud ud via the exchange of two 
gluons.) Therefore, we cannot simply take the NNLO production as calculated and append a Breit-Wigner propagator 
that uses leading or next-to-leading order widths. Instead, we perform the complete calculation at next-to- leading 
order. 

We organize the paper as follows. In Sec. || we calculate the partial widths into leptons or quarks, including the 
effects of the top-quark mass, for the W' boson using arbitrary vector and axial-vector couplings at next-to-leading 
order (NLO) in QCD. We demonstrate that the couplings factorize, and present some numerical results for partial 



widths and branching fractions into a top-quark final state. In Sec. Ill we calculate the fully differential production and 
decay of a W' into any final state at NLO in QCD. We present numerical results for the single-top-quark production 
cross section via s-channel exchange of a W boson. We asses all theoretical uncertainties, and propose a modification 
of the tolerance method for calculating parton distribution function uncertainties. Finally, we present the NLO jet 
distributions, discuss the effect of jet definitions on these distributions, and estimate the limits that may be placed 
on the W mass using data from run I and run II of the Fermilab Tevatron. 

II. NEXT-TO-LEADING ORDER WIDTH 

We divide our evaluation of the W' width into three terms that depend on the final decay products: 

rtot(I^') = T{W' -> tq') + T{W' ^ qq') + T{W' ^ Iv) . (3) 

We separate the partial widths containing a top quark from those containing only massless quarks so that we may 
retain an explicit top-quark mass dependence in the width. 

While the large number of couplings in Eq. (H) appears daunting, a factorization of the couplings occurs in the 
partial widths of Eq. (^ that greatly simplifies the calculation. The leading order partial widths are 

Ti^oiW ^ tq') = -l^\v;^.\\ml,. + mll2) , (4) 

rLo(vr'-99-') = ^l^,VlW, (5) 
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{W ^ W) ^ ^\VL\^'^ , (6) 
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where /3 = 1 — m^/m^,, and we assume = %rn^GF j as in the standard model. Because we have assumed 
there is at most one non-zero mass in the final state, the couplings appear only in the combination 

\9^'hf, I' = |5L sin C y^s, r + W cos C "^^/f/j' • (7) 

Hence, the partial widths of the W boson have the same form as the standard model W boson, with the effect of 
new couplings and GCKM matrix elements absorbed into V^.y .- 

The above factorization of couplings turns out to hold in the next-to-leading order widths as well. The calculations 
of the next-to- leading order partial widths of the W boson were first performed for the massless [Q, one mass 
p3| , and arbitrary mass |]3^ final states many years ago. In order to derive the factorization of arbitrary couplings 
above, however, we rederive the partial widths, but use a newer calculational method whose results are needed in the 



calculation of the fully differential cross section in Sec. Ill 
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We evaluate the next-to-leading order partial widths to quarks by using the phase-space-slicing method with two 
cutoffs [Q. In this method, phase space is divided into three regions: soft (S), hard coUinear (HC), and hard non- 
coUinear (HC). The first two regions are integrated over two-body phase space in n = 4 — 2e dimensions, whereas the 
third region is integrated over three-body phase space in 4 dimensions with cuts. 

The soft region of phase space is defined in the W rest frame by a condition on the energy of the emitted gluon 



Eg<S, 



mw' 



(8) 



where Sg is an arbitrary parameter that must cancel in the final result. The hard region is the complement. Eg > 
dsmw/2. The gluon can only be coUinear with a light quark, and so the coUinear region is defined by comparing 
the invariant mass squared of the gluon and light quark to Scm^,, i.e. the coUinear region is Sqg = 2pq ■ pg < ScTti^,, 
where Pq^g are quark and gluon four- momenta. The final result must be independent of Sc- In practice we retain terms 
logarithmic in Ss or Sc, and drop terms that are linear in the cutoffs. At the end, we take the couplings numerically 
to zero and show the solution contains no residual or Sc dependence. 
The two-body next-to-leading order correction to the width is given by 
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where Cp — 4/3, and Flo is listed in Eqs. (|[ ||). If there is a top quark in the final state, the terms are 
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2\n\Ss)~2\n{S,) 



•-W' 



■In 



'■W 



-2Li2(/3)--ln^ 



1 - In 

2Li2f^ ) -ln"(/3) + 21n(5,,)ln(/3)-ln(5,) (21n(5,) + ^-21n(/3) 



PSc 



Ml 



2 



m 



w 



rat 



- ln(l - A) 



2 + /3 + 2 In 



ml 



-6-ln2(l-A) + 2Li2(A)~2C2 



2(3 - /3) 



(9) 
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where ^2 — 7r^/6, (3 = 1 — ml/rriy^,, and A = 1//3. Note (1 — A) < 0, which means ln(l — A) = ln|l — A|, but 
ln2(l-A) = ln2|l-A|-7r2. 

If both final-state quarks are massless, the formulae simplify to 



Mi 



41n2(,5,) 



= 7-4C2-21n2(,5, 



Ml = -8 -f 4C2 , 



4Li2( J ) -ln(<5e)[41n((5. 



-3], 



V 



0. 



The sum of the terms in the massless case is 

= 2\n^{5s) - ln(5c)[41n(5^ 



3] + 4Li2 



(14) 
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The three-body hard non-coUinear real gluon emission term is evaluated using a Monte Carlo integration in four 
dimensions, subject to the cuts listed above. The integrand is given by 



d(<5r3) = 



1 -, 



(19) 



with 



5]|M3p = -27ra,52c^|y' 



(20) 



and 
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^-r = 2/3m^,(3-/3) 



/ 2 
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+ (3-/3)^ + ^ +2(1-/3), (21) 

where /3 = 1 — rnl/rn^,, s'12 — 2j3i • P2, = 2pi • p3, S23 = 2p2 ■ P3, and pi, p2, and pa are the four-momenta of the 
light quark, top quark, and gluon, respectively. When both quarks are massless /3 = 1, and Eq. ( pl| ) is very simple. 

Once we determine that the factorization of couplings holds for all terms, the partial widths reduce to the expressions 
given in Ref. ||3^] with the replacement of CKM matrix elements by the GCKM matrix elements \ V'\^ defined in Eq. (|^). 
The massless case is simply 

rNL0(M^' - qq') ={l+ ^ VlO m') ■ (22) 



The massive case should reduce to Eq. (12) of Ref. [ p3[ , 

rNLo(W^' - tq') = T^o{W' - tq') + ^^^^^CF\V',\''mw' 



x{/3(l - /3)(/3 - 4) - /32(9 - 5/3) ln(/3) - (1 - /3) (4 + 6/3 - 5/3^) ln(l - /3) 
+/32(3 - /3)[3/2 + 4Li2(l - 13) + 21n(/3) ln(l - (3)]} . (23) 



A. Numerical results 



In order to be of immediate use to experimental analyses at the Fermilab Tevatron, we make definite predictions 
of the W partial and total widths. In all numerical results we use mt = 175 ± 5 GeV, Gp = 1.16639 x 10~^ GeV~^, 
and mw = 80.4 GeV. We use a two-loop running of as as defined in the CTEQ5M1 parton distribution functions 
p6[ . We assume that V^^ is the identity matrix, and for V^'^/ we use the average Cabibbo-Kobayashi-Maskawa (CKM) 
matrix fl^ with the exception that we assume \V^^\ = 1, 

/ 0.9751 0.2215 0.0035 \ 
^CKM = 0.2210 0.9743 0.0410 . (24) 
\ 0.0090 0.0400 1.0000 ) 
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FIG. 1. Dependence ofthe NLO correction to the W width on 5s = 300 x 5^, for mw' = 500 GeV. ST^i^o = 5^2 + SV^,, 
where refers to the HC component, and 5T2 to the sum of virtual, soft and hard-coUinear contributions. The bottom region 
shows an enlargement of the numerical prediction (circles) compared to the analytic prediction (solid line) in Eq. (b3) . 
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Before we discuss our results, we must show that our final result does not depend on the cutoffs we choose. To 
demonstrate this, we fix the ratio of 5s /5c = 300, and plot in Fig. 0the NLO correction to the width of a W' of mass 
500 GeV. In the upper half of the figure we see the logarithmic dependence in 5s and 5c of the two- and three-body 
terms, 5T2 and 5Tj,. In the bottom half of the figure we focus on the sum of the correction terms ^Fnlo for various 
values of 5s (circles), and compare the result to the analytic prediction (solid line) of Eq. (^3|). Once 5s ^ 10"'^ 
{5c ^ 1/3 X 10^^), the result is stable to much less than 1% of the known NLO correction. In practice there is a 
tradeoff between numerical accuracy in canceling the logarithmic divergences of 5s and 5c in the Monte Carlo, and 
the residual power suppressed dependence on these parameters. Hence, we always choose values of 5s and 5c to ensure 
that the residual effects are smaller than the desired Monte Carlo statistical error. 

In Table || we list the LO and NLO partial widths for the decay of the W into light quarks, or tb. Since we assume 
the standard model-like couplings of Eq. (p^), the decays into id and ts are strongly suppressed. However, we remind 
the reader that any model may be restored via the use of Eq. (^ . We present the results for W' masses between 200 
GeV and 1000 GeV, in order to cover the possible reach of experiments at the Tevatron. 



TABLE I. LO and NLO partial widths for W' decays into light quarks, or tb. 



Mass (GeV) 


^qq) (GeV) 


rNLo(W ^ qq') (GeV) 


rLo(W' ^ tb) (GeV) 


rNLo(H/' ^ tb) (GeV) 


200 


1.697 


1.754 


0.129 


0.166 


225 


1.909 


1.973 


0.388 


0.473 


250 


2.121 


2.191 


0.687 


0.811 


275 


2.333 


2.409 


0.993 


1.146 


300 


2.546 


2.628 


1.296 


1.471 


350 


2.970 


3.064 


1.879 


2.084 


400 


3.394 


3.500 


2.432 


2.655 


450 


3.818 


3.936 


2.959 


3.196 


500 


4.243 


4.372 


3.467 


3.715 


550 


4.667 


4.807 


3.961 


4.218 


600 


5.091 


5.243 


4.443 


4.710 


650 


5.515 


5.678 


4.917 


5.192 


700 


5.940 


6.114 


5.384 


5.667 


750 


6.364 


6.549 


5.845 


6.136 


800 


6.788 


6.984 


6.301 


6.600 


850 


7.212 


7.420 


6.754 


7.061 


900 


7.637 


7.855 


7.204 


7.518 


950 


8.061 


8.290 


7.651 


7.972 


1000 


8.485 


8.725 


8.096 


8.425 
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Unless the couplings are suppressed, a left-handed W will decay into either quarks or leptons. A right-handed W, 
however, will decay to leptons only if the mass of the right-handed neutrino (mj^^) is less than niw', or if there is large 
left-right mixing in the neutrino sector. For completeness, we consider both kinematic cases. In Tables || and II] we 
list the LO and NLO W total widths and branching fractions into a tb final state. We assume a top-quark mass of 
175 GcV. The largest uncertainty in the predictions of both the total width and the branching fractions is a result 
of the micertainty in the top-quark mass. Hence, we also show the increase (upper error), or decrease (lower error), 
in the quantities if the top-quark mass is 170 GeV, or 180 GeV, respectively. In Table || we include only the decays 
into quark final states. In Table III we allow for the decay into quarks or leptons as in the standard model. The total 
widths in these Tables are used in Sec. Ill in the NLO calculation of the W contribution to the single-top-quark cross 
section. 



TABLE II. Total width of the W', and its branching ratio into tb, at LO and NLO in QCD, when the decay to leptons is 
not allowed. Errors are due to the top-quark mass uncertainty, mt = 175 =F 5 GeV. 



Mass (GeV) 



Tlo (GeV) 



BRlo(IF' -> fb) 



Lnlo (GeV) 



BRNLo(M^'->tb) 



200 


3.523 


+0 
-0 


049 
043 


0.0366 


+0.0132 
-0.0119 


3.675 


+0.060 
-0.053 


0.0452 


+0.0152 
-0.0140 


225 


4.206 


+0 
-0 


064 
061 


0.0923 


+0.0136 
-0.0135 


4.419 


+0.073 
-0.071 


0.1071 


+0.0144 
-0.0146 


250 


4.930 


+0 
-0 


068 
067 


0.1394 


+0.0117 
-0.0119 


5.193 


+0.074 
-0.074 


0.1561 


+0.0119 
-0.0122 


275 


5.660 


+0 
-0 


067 
068 


0.1755 


+0.0097 
-0.0100 


5.965 


+0.071 
-0.072 


0.1922 


+0.0095 
-0.0099 


300 


6.388 


+0 
-0 


065 
066 


0.2030 


+0.0081 
-0.0083 


6.726 


+0.067 
-0.068 


0.2187 


+0.0077 
-0.0080 


350 


7.819 


+0 
-0 


059 
060 


0.2404 


+0.0057 
-0.0059 


8.211 


+0.058 
-0.060 


0.2537 


+0.0053 
-0.0055 


400 


9.220 


+0 
-0 


053 
054 


0.2637 


+0.0042 
-0.0044 


9.654 


+0.051 
-0.052 


0.2750 


+0.0038 
-0.0040 


450 


10.595 


+0 
-0 


048 
049 


0.2792 


+0.0032 
-0.0033 


11.068 


+0.045 
-0.046 


0.2888 


+0.0029 
-0.0030 


500 


11.952 


+0 
-0 


043 
044 


0.2901 


+0.0026 
-0.0027 


12.458 


+0.040 
-0.041 


0.2982 


+0.0022 
-0.0023 


550 


13.294 


+0 
-0 


040 

041 


0.2979 


+0.0021 
-0.0022 


13.832 


+0.036 
-0.037 


0.3049 


+0.0018 
-0.0019 


600 


14.625 


+0 
-0 


036 
037 


0.3038 


+0.0017 
-0.0018 


15.195 


+0.032 
-0.034 


0.3099 


+0.0015 
-0.0015 


650 


15.947 


+0 

-0 


034 
035 


0.3083 


+0.0015 
-0.0015 


16.548 


+0.030 
-0.031 


0.3137 


+0.0012 
-0.0013 


700 


17.263 


+0 
~0 


031 
032 


0.3119 


+0.0012 
-0.0013 


17.894 


+0.027 
-0.028 


0.3167 


+0.0010 
-0.0011 


750 


18.572 


+0 
-0 


029 
030 


0.3147 


+0.0011 
-0.0011 


19.234 


+0.025 
-0.026 


0.3190 


+0.0009 
-0.0009 


800 


19.878 


+0 
-0 


027 
028 


0.3170 


+0.0009 
-0.0010 


20.569 


+0.023 
-0.024 


0.3209 


+0.0008 
-0.0008 


850 


21.179 


+0 
-0 


026 

026 


0.3189 


+0.0008 
-0.0009 


21.900 


+0.021 
-0.023 


0.3224 


+0.0007 
-0.0007 


900 


22.477 


+0 
-0 


024 
025 


0.3205 


+0.0007 
-0.0008 


23.227 


+0.021 
-0.021 


0.3237 


+0.0006 
-0.0006 


950 


23.773 


+0 
-0 


023 
024 


0.3218 


+0.0007 
-0.0007 


24.552 


+0.020 
-0.019 


0.3247 


+0.0005 
-0.0005 


1000 


25.066 


+0 
~0 


022 
023 


0.3230 


+0.0006 
-0.0006 


25.875 


+0.018 
-0.018 


0.3256 


+0.0005 
-0.0005 
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TABLE III. Total width of the W, and its branching ratio into tb, at LO and NLO in QCD, when decays to quarks or 
leptons are both included. Errors are due to the top-quark mass uncertainty, mt = 175 =F 5 GeV. 



Mass (GeV) Flo (GeV) BRj^p {W ^ tb) Fnlo (GeV) BRnlo(W^' ^ tb) 



200 


5.220 


+0.049 
-0.043 


0.0247 


+0.0091 
-0.0081 


5.372 


+0.060 
-0.053 


0.0309 


+0.0106 
-0.0097 


225 


6.116 


+0.064 
-0.061 


0.0635 


+0.0097 
-0.0095 


6.328 


+0.073 
-0.071 


0.0748 


+0.0105 
-0.0105 


250 


7.051 


+0.068 
-0.067 


0.0974 


+0.0086 
-0.0087 


7.314 


+0.074 
-0.074 


0.1109 


+0.0089 
-0.0091 


275 


7.994 


+0.067 
-0.068 


0.1243 


+0.0073 
-0.0075 


8.298 


+0.071 
-0.072 


0.1381 


+0.0073 
-0.0076 


300 


8.933 


+0.065 
-0.066 


0.1451 


+0.0062 
-0.0064 


9.272 


+0.067 
-0.068 


0.1587 


+0.0060 
-0.0063 


350 


10.789 


+0.059 
-0.060 


0.1742 


+0.0046 
-0.0046 


11.181 


+0.058 
-0.060 


0.1863 


+0.0042 
-0.0044 


400 


12.614 


+ 0.053 
-0.054 


0.1928 


+0.0034 
-0.0035 


13.049 


+0.051 
-0.052 


0.2035 


+0.0031 
-0.0032 


450 


14.414 


+0.048 
-0.049 


0.2053 


+0.0026 
-0.0027 


14.886 


+0.045 
-0.046 


0.2147 


+0.0023 
-0.0024 


500 


16.195 


+0.043 
-0.044 


0.2141 


+0.0021 
-0.0022 


16.701 


+0.040 
-0.041 


0.2224 


+0.0019 
-0.0019 


550 


17.961 


+0.040 
-0.041 


0.2205 


+0.0017 
-0.0018 


18.499 


+0.036 
-0.037 


0.2280 


+0.0016 
-0.0015 


600 


19.716 


+0.036 
-0.037 


0.2253 


+0.0014 
-0.0015 


20.287 


+0.032 
-0.034 


0.2322 


+0.0012 
-0.0013 


650 


21.463 


+ 0.034 
-0.035 


0.2291 


+ 0.0012 
-0.0012 


22.064 


+ 0.030 
-0.031 


0.2353 


+ 0.0010 
-0.0011 


700 


23.202 


+0.031 
-0.032 


0.2320 


+0.0010 
-0.0011 


23.834 


+0.027 
-0.028 


0.2378 


+0.0009 
-0.0009 


750 


24.936 


+0.029 
-0.030 


0.2344 


+0.0009 
-0.0009 


25.598 


+0.025 
-0.026 


0.2397 


+0.0007 
-0.0008 


800 


26.666 


+0.027 
-0.028 


0.2363 


+0.0008 
-0.0008 


27.357 


+0.023 
-0.024 


0.2413 


+0.0006 
-0.0007 


850 


28.391 


+0.026 
-0.026 


0.2379 


+0.0007 
-0.0007 


29.113 


+0.021 
-0.023 


0.2425 


+0.0006 
-0.0006 


900 


30.114 


+0.024 
-0.025 


0.2392 


+0.0006 
-0.0006 


30.864 


+0.021 
-0.021 


0.2436 


+0.0005 
-0.0005 


950 


31.834 


+0.023 
-0.024 


0.2403 


+0.0006 
-0.0006 


32.613 


+0.020 
-0.019 


0.2444 


+0.0005 
-0.0004 


1000 


33.551 


+0.022 
-0.023 


0.2413 


+0.0005 
-0.0005 


34.360 


+0.018 
-0.018 


0.2452 


+0.0004 
-0.0004 
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It is not surprising that top-quark threshold effects cause a large uncertainty in the branching fraction for W masses 
less than 300 GeV. This also appears in the large increase in branching fraction at NLO over the LO branching 
fraction. With 33% changes in both the mass variation and NLO correction, perturbation theory is somewhat suspect 
if m\Y' ^ 200 GeV. However, the effect is less than 10% by 225 GeV, and rapidly vanishes as m^, ^ mf. The 
branching fraction into a top-quark final state is nearly saturated at 1/3 (1/4) when myy' ~ 500 GeV for the quark- 
only ( quark-plus- I cpton) model of decays. Hence, a large fraction of W events should produce a top-quark in the 



final state. In Sec. HI A we see this effect on the single-top-quark cross section. 



III. NEXT-TO-LEADING ORDER CROSS SECTION 

The analytic form of the differential production and decay of a W at next-to-leading order is very similar to that 
of s-channel single-top-quark production. The complete calculation of differential single-top-quark production using 
the phase space slicing method appears in Ref. . We generalize the calculation in that paper to the production of 
one massive particle through a W boson (such that single-top-quark production is a special case) . We follow closely 
the notation of Ref. [Q, but retain arbitrary couplings in the vertices. We note that the <-channel exchange of a 
virtual W is suppressed by at least l/mfy,, and hence do not consider it here. However, the analytic expressions and 
factorizations that follow are valid for this channel as well. The analytic formulae for the i-channel exchange may be 
obtained by simple crossing. 

At leading order, the production of a heavy quark (lepton) may be written schematically as 

u(pi)J(p2) ^ W'+ ^ bip^KpiMpsHpi)] , (25) 

where pi + P2 = Ps + Pi, ud represents all possible parton fluxes, and the heavy particle in the final state (t or lyji) 
has four-momentum p^ and mass m. 

The leading-order fully differential spin-averaged partonic cross section can be written as 

dao - ^^|A^o|'dPS2, (26) 
where s — (pi + P2)'^ is the partonic center-of-momentum energy squared, and the two body phase space is given by 

The matrix element squared, summed (averaged) over final- (initial-) state spin and color states, for production of 
one massive particle is 

2^\Mo\ -— —2 — — , (28) 

where t — {pi — pa)^, u — {p2 — Ps)^- The number of colors Nc = 3, and m — mt, for a final state with a top quark. 
For the production of a massive neutrino m — m^j^ and Nc = 1. If the final state is massless, simply set m = and 
use Nc = 3 (1) for quarks (leptons). The functions Rt and are given by 

Rt = {\V,f + \Af){\Vff + \Aff)+4Re{V.AnRe{VfA/) , (29) 
Ru - {m' + M{\Vff + \Aff)-4Re{V,AnMVfA/), (30) 

where Vij and Aij are the vector and axial-vector couplings for the initial- and final-state vertices, respectively. In 
the notation of Eq. (||), 

Iffi? cos C . I ^ + sin C l^/'f . I ^ 
\V\^ + \A\^ = \^ ^ , (31) 

2Re{VA*) = ^— ' ^ ' ^_JiIiL . (32) 

For right- or left-handed W bosons \V\ = i?„ = 0, Rt reduces to Rt = (.g^l^/Pl V^P)/8, and we recover the 
standard model s-channel single-top-quark cross section published in Refs. | p7| , p8[ , including the Breit-Wigner term, 
up to the GCKM couphngs. 
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We calculate the fully differential next-to-leading order cross section with the same method used in Sec. ||. Again, 
the phase space is divided into a hard and soft region, where the soft region is now defined in the partonic center-of- 
momentum frame by a condition of the emitted gluon energy: 



Eg < 5s ^ , (33) 



where Sij = [pi + pj)'^. The hard region is the complement, Eg > Ss ^/syill. A region is coUinear if an invariant in 
a propagator is less than (5cSi2, i.e. the region is coUinear if in a denominator we obtain = (pi < (5cSi2, or 

= [Pi -PoY < 5cSl2- 

At next-to-leading order, color conservation forbids the exchange of a single gluon between the initial and final 
states. This is convenient when organizing the solution, because the leptonic cross section only has QCD corrections 
in the initial state. Hence, we list separately the initial and final state corrections. 

The two-body NLO correction to the partonic cross section contains 

dcT2 ^^CFml + Ml -f Ml)da^ + dav], (34) 

where Cp = 4/3, and day is the part of the correction to the vertex containing the top quark that is not proportional 
to the LO cross section. If the final state does not contain a top quark then day — 0; otherwise, 

d^v = {Rt^'Ru)-- f-^ ^-—dVS^. (35) 

The initial-state corrections appearing in Eq. (^) are 

Ml, = 4ln'{Ss), (36) 
M2^ = [41n(5,) + 3]lnf4V (37) 



M^, = -8 + 4C2 , (38) 

where C2 = 7r^/6, and p,F is the factorization scale. If there is a top quark in the final state, Mgj:, Mqj, and Myj: 
are the same as those listed in Eqs. (p^|-|T2|), with the replacements m^, ^ s, /3 ^ 1 — rri^/s, and X = 1/(3. If both 
quarks in the final state are massless, then Mj is given by Eq. (18). 
The complete two-body cross section is given by 

- E/ dx,dx2{f^'{xi,fiF)f,'''{x2,fiF)d<jf + [/f^(xi,/i;^)/f^(x2,/iF) + f!^H^uf^F)f"'{x2,f^F)]d<jS'} , (39) 

where a, b sum over all quark- antiquark luminosities, f^p is the factorization scale. Hi 2 are the initial-state hadrons, 
and the f{x, fi) are NLO parton distribution functions. The / functions are introduced to compensate for a difference 
between the limits of integration used in the phase space slicing calculation of the initial-state coUinear singularities 
and the modified minimal subtraction MS scheme used in the NLO parton distribution functions. The / functions 
are given in Ref. Q . 

The three-body hard non-coUinear terms are evaluated using a Monte Carlo integration in four dimensions, subject 
to the cuts listed above. The cross section is given by 



Cr3=-V / dxidX2^^^ [ _{f^J{xi,^iF)f"Hx2,^^F)'i'qq 



a,b " " -^HC 

+ [/,^i(xi,/iF)/fMa^2,MF) +/fn2;i,MF)/gT'(2^2,AiF)]*gs}dPS3. (40) 

If we label the momenta in the three-body processes 

■^gq : u{pi)d{p2) -> b{p3)t{p4)9{P5) , (41) 
'fgg : u{pi)g{p2) b{p3)t{p4)d{p5) , (42) 

then functions ^^g and "^gg are given by 
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2Cf 



Rt - 



^24 (^13 + ^35) ^13^14 ^13 (^24 + ■S45) ^23^24 



■^12 [^13 (2^24 + ^45) + '^35^24] 



^25 
Ru iPl 



-P2] 



2Cf 



2 

w'^ w 



Rt - 



^24 (^13 + ^15) ~ ^13^23 *13 (^24 + ^25) ~ /S45) — 



'45 



S34 [^13 (2*24 + ^25) + ^15*24] 



'35 ''45 



Ru (pi 



P2] 



Rt 



■S45 (*13 + ^23) ~ ^13*14 , *13 (*45 + ^24) ~ ■^35*45 



''25 



'Ru 



''23 ^ 

■^12 

*15 [^45 (2*13 + ^23) + ^13*24] 
^12*25 

'35 (^'14 + ^24) ~ *13*i4 I ^'14 (^35+^23) 



'35''45 



"25 



*15 [■S35 (2*14 + *24) + *'l4*23] 



S]^2* 



25 



(43) 



(44) 



where s', — 



ml, — tij — m^, and all other terms may be obtained by crossing. Final-state QCD corrections 



are limited to the second term of '^qq, and thus this term does not appear in corrections to the leptonic final state. If 
the final state is massless, then the solution is recovered by setting rrit = 0, and noting that all primed invariants are 
equal to their unprimed counterparts. 



A. Single-top-quark production via ^ 



In Sec. II A we saw that W' bosons tend to have a large branching fraction into top quarks. This observation 



leads us to consider the effect of such a W on the size of the single-top-quark cross section at hadron colliders. In 
particular, we show the cross sections are large enough to improve the mass limits on W' bosons using data from run 
I of the Fermilab Tevatron. The W boson affects the single-top-quark cross section through the three channels shown 
in Fig. ||. We concentrate on the s-channel production, because of the enhancement from the W resonance. We do 
not present numerical results for the i-channel or associated production of a W boson, because the cross sections for 
these channels are negligible for the masses we consider. 




(a) (b) (c) 

FIG. 2. Representative Feynman diagrams for single-top-quark production involving a W' boson: (a) s-channel production, 
(b) i-channel production, and (c) W' — t associated production. 



In order to make definitive predictions, we stick to the assumption that the W has purely right- or left-handed 
interactions. We choose the standard model CKM matrix as in Eq. (|2j), with the und erstanding that the final results 
may be scaled using the factorization of couplings in the last section. As in Sec. [I A , we choose m* = 175 ± 5 GeV, 
Gf = 1.16639 X 10~^ GeV~^, and mw = 80.4 GeV. For leading-order cross sections we use CTEQ5L leading-order 
parton distribution functions (PDFs). We use a two- loop running of a^, and CTEQ5M1 PDFs for NLO cross sections. 
We use the LO and NLO widths calculated in Sec. II A in the LO and NLO cross sections, respectively, so that all 



terms are calculated at the same order and with the same method. We set both factorization and renormalization 
scales equal to mw' , since most of the cross section will be produced near resonance. 

We present the leading- and next-to-leading order cross sections for single-top-quark production, via ^ bosons 

of various masses, at run I of the Tevatron (a pp collider with y/S = 1.8 TeV) in Tables and 0. In Table nM we 
perform the calculation under the assumption that vji is too heavy for the W' to decay to leptons, and thus use the 
widths of Tabic 0. In Table we assume all decays are allowed, and use the widths of Table |l|. The cross sections 
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for just top-quark or antitop-quark production are one-half of the hsted results. The corresponding si ngle- top-quark 
cross sections a t a r un II of the Tevatron (a pp colhder with ^/S = 1.96 TeV) appear in Tables VI and VI] . 

In Tables |v|- |VIl| we show the uncertainties in the cross sections due to the variation of the scale between mw' /2 
and 2mw' , the variation of the top-quark mass over mt = 175 ± 5 GeV including the effects of the change in width 
listed in Tables ^ and [II, and uncertainties from the parton distribution functions (using a new method described 
in Sec. Ill B ). Additional uncertainties come from the use of a vanishing bottom-quark mass, which overestimates 



the cross section by +1.5% at 200 GeV, but only -1-0.4% by 225 GeV, and is a negligible effect at higher W' masses. 
Extrapolating from the results of Ref. we estimate that Yukawa corrections cause a shift of less than —1% in 
the cross section. The uncertainty from Monte Carlo statistics is 0.03%. There is a ±4% uncertainty in the NLO 
correction term from Ug = 0.118 ± 0.005. All uncertainties listed above are added in quadrature, and are presented 
for easy reference in the last column of each section in Tables [V-VII as a percentage error. If we had used the 



narrow width approximation, and multiplied the NLO W cross section times NLO branching fraction to tb instead 
of evaluating the full matrix element, there would be an additional ±1-2% systematic uncertainty. 



TABLE IV. LO and NLO cross sections in (pb) for pp W^^ ^ tb + tb at run I of the Tevatron, ^/5 = 1.8 TeV, when 
the decay to leptons is not allowed. Scale, top-quark mass, and PDF uncertainties are also listed in (pb). The last column in 
each section lists the total theoretical uncertainty in (%), where all uncertainties in the text and in this Table are added in 
quadrature. 



iViaSS l^VjrG V ) 


o"LO (pb) 


Sai^o{fJ-,Smt, PDF) (pb) 




0"NLO (pb) 


SaKi.o{^,Smt, PDF) (pb) 5a'^°j^o (%) 


200 


32.24 


+2.0 
-1.8 


+ 10.9 
-9.8 


+ 1.3 
-1.2 


+34.7 
-31.2 


50.90 


+3.0 
-2.4 


+16 
-15 


+2.1 
-1.9 


+32.3 
-30.2 


225 


51.45 


+3.6 

-3.2 


+7.5 
-7.4 


+2.3 
-2.0 


+ 16.8 
-16.2 


77.24 


+4.3 
-3.6 


+10 
-10 


+3.4 
-3.0 


+ 14.9 
-14.4 


250 


52.29 


+4.1 

-3.6 


+4.4 
-4.5 


+2.5 
-2.1 


+ 12.5 
-11.8 


76.10 


+4.1 
-3.6 


+5.9 
-6.1 


+3.6 
-3.0 


+ 10.7 
-10.2 


275 


45.58 


+3.9 
-3.4 


+2.5 
-2.6 


+2.3 
-1.9 


+ 11.4 
-10.3 


64.92 


+3.4 
-3.1 


+3.3 
-3.4 


+3.3 
-2.7 


+9.0 
-8.4 


300 


37.18 


+3.4 
-2.9 


+ 1.5 
-1.6 


+2.0 
-1.6 


+ 11.4 
-9.9 


52.11 


+2.8 
-2.6 


+ 1.9 
-2.0 


+2.8 
-2.2 


+8.6 
-7.7 


350 


22.82 


+2.4 
-2.0 


+0.56 
-0.58 


+ 1.4 
-1.1 


+ 12.5 
-10.4 


31.20 


+ 1.7 
-1.6 


+0.70 
-0.73 


+ 1.9 
-1.4 


+8.6 
-7.4 


400 


13.46 


+1.5 
-1.3 


+0.23 
-0.23 


+0.91 
-0.68 


+ 13.2 
-11.1 


18.04 


+ 1.03 
-0.98 


+0.28 
-0.29 


+ 1.22 
-0.91 


+9.1 
-7.7 


450 


7.836 


+0.98 
-0.82 


+0.10 
-0.10 


+0.59 
-0.44 


+ 14.7 
-12.0 


10.31 


+0.62 
-0.59 


+0.12 
-0.12 


+0.78 
-0.58 


+9.8 
-8.2 


500 


4.537 


+0.61 
-0.51 


+0.045 
-0.046 


+0.39 
-0.28 


+ 16.0 
-12.9 


5.873 


+0.38 
-0.36 


+0.054 
-0.056 


+0.50 
-0.37 


+ 10.8 
-9.0 


550 


2.617 


+0.38 
-0.31 


+0.022 
-0.022 


+0.26 
-0.18 


+ 17.6 
-13.8 


3.335 


+0.23 
-0.21 


+0.026 
-0.027 


+0.33 
-0.23 


+ 12.2 
-9.5 


600 


1.503 


+0.23 
-0.19 


+0.011 
-0.011 


+0.17 
-0.12 


+ 19.1 
-15.0 


1.889 


+0.14 
-0.13 


+0.013 
-0.013 


+0.21 
-0.15 


+ 13.4 
-10.6 


650 


0.859 


+0.14 
-0.11 


+0.006 
-0.006 


+0.11 
-0.07 


+20.8 
-15.2 


1.066 


+0.08 
-0.08 


+0.007 
-0.007 


+0.13 
-0.09 


+ 14.4 

-11.4 
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TABLE V. LO and NLO cross sections in (pb) for pp W^^ ^ ^ tb + tb at run I of the Tevatron, \/S = 1.8 TeV, when 
decays to quarks or leptons are both included. Scale, top-quark mass, and PDF uncertainties are also listed in (pb). The last 
column in each section lists the total theoretical uncertainty in (%), where all uncertainties in the text and in this Table are 
added in quadrature. 



Mass (GeV) 


o-LO (pb) 


Jctlo (At,5mt, PDF) (pb) 


Sal^' (%) 


cnlo (pb) 


5cTNLo(/i,5mt, PDF) (pb) Sa^to (%) 


200 


21.91 


+ 1.4 
-1.3 


+7.3 
-6.5 


+0.91 
-0.84 


+34.2 
-30.5 


34.94 


+ 2.0 
-1.7 


+ 11.1 
-10.1 


+ 1.5 
-1.3 


+32.7 
-29.6 


225 


35.05 


+2.5 
-2.2 


+5.2 

-5.1 


+1.5 

-1.4 


+17.0 
-16.4 


53.36 


+2.9 

-2.5 


+7.5 

-7.4 


+2.4 

-2.1 


+15.8 
-15.3 


250 


36.15 


+2.8 
-2.6 


+3.2 
-3.2 


+ 1.7 
-1.4 


+12.7 
-11.9 


53.40 


+2.9 
-2.5 


+4.4 
-4.6 


+2.5 
-2.1 


+ 11.1 
-10.6 


275 


31.92 


+2.7 
-2.4 


+ 1.9 
-1.9 


+1.6 
-1.3 


+11.6 
-10.5 


46.14 


+2.4 
-2.2 


+2.5 
-2.6 


+2.3 
-1.9 


+9.2 
-8.6 


300 


26.31 


+2.4 
-2.1 


+ 1.2 
-1.2 


+1.4 
-1.1 


+11.6 
-10.1 


37.41 


+2.0 
-1.8 


+1.5 
-1.6 


+2.0 
-1.6 


+8.7 
-7.9 


350 


16.40 


+ 1.7 
-1.4 


+t).44 
-0.46 


+0.98 
-0.75 


+ 12.3 
-10.1 


22.72 


+ 1.2 
-1.2 


+0.56 
-0.58 


+ 1.4 
-1.0 


+8.6 
-7.5 


400 


9.773 


+ 1.12 
-0.95 


+0.18 
-0.19 


+0.66 
-0.49 


+ 13.5 
-11.2 


13.26 


+0.75 
-0.72 


+0.23 
-0.24 


+0.89 
-0.67 


+9.1 
-7.8 


450 


5.735 


+0.72 
-0.60 


+0.081 
-0.083 


+0.44 
-0.32 


+14.8 
-12.0 


7.641 


+0.46 
-0.44 


+0.10 
-0.10 


+0.58 
-0.43 


+9.9 
-8.3 


500 


3.343 


+0.45 
-0.37 


+0.038 
-0.039 


+0.29 
-0.21 


+16.1 
-12.8 


4.380 


+0.28 
-0.27 


+0.046 
-0.047 


+0.38 
-0.27 


+10.9 
-8.9 


550 


1.941 


+0.28 
-0.23 


+ 0.018 
-0.019 


+ 0.19 
-0.14 


+ 17.5 
-13.9 


2.503 


+0.17 
-0.16 


+0.022 
-0.023 


+ 0.24 
-0.17 


+ 11.9 
-9.5 


600 


1.122 


+0.17 

-0.14 


+0.009 

-0.009 


+0.12 
-0.09 


+18.6 
-14.9 


1.427 


+0.10 
-0.10 


+0.011 
-0.011 


+0.16 
-0.11 


+13.3 
-10.5 


650 


0.646 


+0.11 
-0.08 


+0.005 
-0.005 


+0.08 
-0.06 


+21.1 
-15.5 


0.812 


+0.06 
-0.06 


+0.006 
-0.006 


+0.10 
-0.07 


+14.4 
-11.5 
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TABLE VI. LO and NLO cross sections in (pb) for pp -» W^^ ^ -» tfe + i6 at run II of the Tevatron, -/S = 1.96 TeV, when 
the decay to leptons is not allowed. Scale, top-quark mass, and PDF uncertainties are also listed in (pb). The last column in 
each section lists the total theoretical uncertainty in (%), where all uncertainties in the text and in this Table are added in 
quadrature. 



Mass (GeV) 


o-LO (pb) 


Saj^oilijmt, PDF) (pb) 


StI^' (%) 


fNLO (pb) 


<5iTNLo(At,5mt, PDF) (pb) 




200 


36.39 


+2.1 
-1.9 


+ 12.3 
-11.0 


+ 1.5 
-1.4 


+34.6 
-31.0 


57.39 


+3.2 
-2.6 


+ 18.4 
-16.8 


+2.3 
-2.2 


+33.1 
-30.2 


225 


58.56 


+3.8 
-3.4 


+8.5 

-8.4 


+2.5 

-2.2 


+ 16.6 

-15.9 


87.86 


+4.6 
-3.9 


+11.8 
-11.9 


+3.7 
-3.4 


+15.6 

-15.3 


250 


60.13 


+4.3 
-3.8 


+5.0 
-5.1 


+2.7 
-2.3 


+ 11.9 
-11.2 


87.52 


+4.6 
-3.9 


+6.8 
-6.9 


+3.9 
-3.4 


+ 11.1 
-10.7 


275 


53.00 


+4.1 
-3.6 


+3.0 
-3.0 


+2.5 
-2.1 


+ 10.7 
-9.7 


75.57 


+3.8 
-3.4 


+3.9 
-4.0 


+3.6 
-3.0 


+9.6 
-9.0 


300 


43.78 


+3.7 
-3.2 


+1.8 
-1.8 


+2.2 
-1.8 


+ 10.7 
-9.3 


61.46 


+3.1 
-2.8 


+2.3 
-2.4 


+3.1 
-2.5 


+9.0 
-8.3 


350 


27.63 


+2.6 
-2.3 


+0.68 

-0.70 


+ 1.6 
-1.2 


+ 11.3 
-9.7 


37.90 


+ 1.9 
-1.8 


+0.85 
-0.88 


+2.1 
-1.7 


+8.8 
-8.1 


400 


16.81 


+1.8 
-1.6 


+0.28 
-0.29 


+ 1.05 
-0.80 


+ 12.5 
-10.3 


22.65 


+ 1.2 
-1.1 


+0.35 
-0.36 


+ 1.4 
-1.1 


+9.3 
-8.2 


450 


10.13 


+1.2 
-1.0 


+0.13 
-0.13 


+0.70 
-0.62 


+ 13.8 
-11.2 


13.43 


+0.76 
-0.72 


+0.16 
-0.16 


+0.93 
-0.70 


+9.9 
-8.6 


500 


6.092 


+0.76 
-0.64 


+0.060 
-0.061 


+0.47 
-0.36 


+ 14.7 
-12.0 


7.953 


+0.47 
-0.45 


+0.072 
-0.074 


+0.62 
-0.45 


+10.7 
-9.0 


550 


3.661 


+0.49 
-0.41 


+0.030 
-0.030 


+0.32 
-0.23 


+ 16.0 
-12.9 


4.710 


+0.30 
-0.28 


+ 0.035 
-0.036 


+0.41 
-0.30 


+11.6 
-9.7 


600 


2.197 


+0.32 
-0.26 


+0.015 
-0.016 


+0.21 

-0.15 


+ 17.4 
-13.7 


2.790 


+0.18 
-0.18 


+0.018 

-0.019 


+0.27 
-0.19 


+ 12.4 

-10.3 


650 


1.316 


+0.20 
-0.16 


+0.0081 
-0.0082 


+0.14 
-0.10 


+ 18.6 
-14.4 


1.650 


+0.11 
-0.11 


+0.010 
-0.010 


+0.18 
-0.13 


+ 13.4 
-11.1 


700 


0.7855 


+0.13 
-0.10 


+0.0045 
-0.0046 


+0.096 
-0.067 


+20.6 
-16.3 


0.974 


+0.071 
-0.069 


+0.0064 
-0.0055 


+0.12 
-0.08 


+14.9 
-11.6 


750 


0.4671 


+0.080 
-0.063 


+0.0026 
-0.0026 


+0.063 
-0.044 


+21.8 
-16.5 


0.573 


+0.043 
-0.043 


+0.0032 
-0.0031 


+0.077 
-0.054 


+15.9 
-12.7 


800 


0.2767 


+0.049 
-0.039 


+0.0015 
-0.0015 


+0.041 
-0.028 


+23.1 
-17.4 


0.337 


+0.027 
-0.027 


+0.0019 
-0.0019 


+0.050 
-0.034 


+ 17.4 
-13.5 


850 


0.1634 


+0.030 
-0.024 


+0.0010 
-0.0010 


+0.026 
-0.018 


+24.3 
-18.4 


0.198 


+0.016 
-0.016 


+0.0012 
-0.0012 


+0.032 
-0.022 


+ 18.5 
-14.4 


900 


0.0963 


+0.018 
-0.014 


+0.0006 
-0.0006 


+0.017 
-0.011 


+25.7 
-18.5 


0.116 


+0.010 
-0.010 


+0.0008 
-0.0008 


+0.020 
-0.014 


+19.7 
-16.4 


950 


0.0569 


+0.011 
-0.009 


+0.0004 
-0.0004 


+0.010 
-0.007 


+26.1 
-20.1 


0.0688 


+0.0069 
-0.0060 


+0.0006 
-0.0006 


+0.0127 
-0.0084 


+20.8 
-16.6 


1000 


0.0339 


+0.006 
-0.005 


+0.0003 
-0.0003 


+0.0064 
-0.0042 


+25.9 
-19.3 


0.0413 


+0.0036 
-0.0036 


+0.0004 
-0.0004 


+0.0079 
-0.0051 


+21.4 
-15.7 



(%) 
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TABLE VII. LO and NLO cross sections in (pb) for pp W^^ L^tb + tb at run II of the Tevatron, = 1.96 TeV, when 
decays to quarks or leptons are both included. Scale, top-quark mass, and PDF uncertainties are also listed in (pb). The last 
column in each section lists the total theoretical uncertainty in (%), where all uncertainties in the text and in this Table are 
added in quadrature. 



Mass (GeV) 


o-LO (pb) 


Saj^oilijmt, PDF) (pb) 


StI^' (%) 


fNLO (pb) 


<5iTNLo(At,5mt, PDF) (pb) 




200 


24.75 


+1.4 
-1.3 


+8.3 
-7.3 


+0.99 
-0.94 


+34.6 
-31.0 


39.43 


+2.2 
-1.8 


+ 12.5 
-11.3 


+ 1.6 
-1.6 


+32.7 
-29.6 


225 


39.92 


+2.6 

-2.3 


+6.0 

-5.8 


+1.7 

-1.5 


+ 16.6 

-15.9 


60.74 


+3.2 
-2.7 


+8.5 
-8.4 


+2.6 
-2.3 


+16.1 

-15.6 


250 


41.58 


+3.0 
-2.6 


+3.7 
-3.7 


+ 1.9 
-1.6 


+ 11.9 
-11.2 


61.43 


+3.1 
-2.7 


+5.0 
-6.1 


+2.7 
-2.4 


+ 11.3 
-11.0 


275 


37.13 


+2.9 
-2.6 


+2.2 
-2.3 


+1.8 
-1.6 


+ 10.7 
-9.7 


53.73 


+2.7 
-2.4 


+3.0 
-3.0 


+2.6 
-2.1 


+9.8 
-9.1 


300 


30.99 


+2.6 
-2.3 


+1.4 
-1.4 


+1.6 
-1.3 


+ 10.7 
-9.3 


44.14 


+2.2 
-2.0 


+1.8 
-1.8 


+2.2 
-1.8 


+9.1 
-8.4 


350 


19.85 


+ 1.9 
-1.6 


+ 0.54 
-0.55 


+ 1.11 
-0.87 


+ 11.3 
-9.7 


27.61 


+ 1.4 
-1.3 


+0.68 
-0.70 


+ 1.5 
-1.2 


+8.8 

-8.0 


400 


12.20 


+ 1.3 
-1.1 


+0.23 
-0.23 


+0.76 
-0.68 


+ 12.6 
-10.3 


16.65 


+0.88 
-0.84 


+0.28 
-0.29 


+ 1.04 
-0.79 


+9.3 
-8.2 


450 


7.411 


+0.86 
-0.73 


+0.10 
-0.11 


+0.52 
-0.38 


+ 13.8 
-11.2 


9.941 


+0.66 
-0.53 


+0.13 
-0.13 


+0.69 
-0.51 


+9.9 
-8.6 


500 


4.484 


+0.56 
-0.47 


+0.050 
-0.061 


+0.35 
-0.26 


+ 14.7 
-12.0 


5.922 


+0.36 
-0.33 


+0.061 
-0.062 


+0.46 
-0.34 


+10.6 
-9.1 


550 


2.709 


+0.36 
-0.30 


+0.025 
-0.026 


+0.24 
-0.17 


+ 16.0 
-12.9 


3.526 


+ 0.22 
-0.21 


+ 0.030 
-0.031 


+0.31 
-0.22 


+11.6 
-9.6 


600 


1.635 


+0.23 
-0.19 


+0.013 
-0.013 


+0.16 
-0.11 


+ 17.4 
-13.7 


2.100 


+0.14 
-0.13 


+0.016 
-0.016 


+0.20 
-0.15 


+ 12.4 
-10.3 


650 


0.9849 


+0.15 
-0.12 


+0.0071 
-0.0072 


+0.108 
-0.076 


+ 18.6 
-14.4 


1.250 


+0.086 
-0.084 


+0.0086 
-0.0087 


+0.14 
-0.10 


+ 13.8 
-11.3 


700 


0.5919 


+0.095 
-0.076 


+0.0040 
-0.0040 


+0.072 
-0.060 


+20.6 
-16.3 


0.7428 


+0.064 
-0.053 


+0.0048 
-0.0049 


+0.090 
-0.063 


+14.7 
-11.8 


750 


0.3547 


+0.060 
-0.047 


+0.0023 
-0.0023 


+0.048 
-0.033 


+21.8 
-16.5 


0.4410 


+0.033 
-0.033 


+0.0029 
-0.0029 


+0.069 
-0.041 


+15.9 
-12.6 


800 


0.2122 


+0.037 
-0.029 


+0.0014 
-0.0014 


+0.031 
-0.022 


+23.1 
-17.4 


0.2619 


+0.020 
-0.020 


+0.0018 
-0.0018 


+0.039 
-0.027 


+ 17.2 
-13.5 


850 


0.1268 


+0.023 
-0.018 


+0.0009 
-0.0009 


+0.020 
-0.014 


+24.3 
-18.4 


0.1558 


+0.013 
-0.013 


+0.0012 
-0.0012 


+0.026 
-0.017 


+ 18.6 
-14.4 


900 


0.0759 


+0.014 
-0.011 


+0.0006 
-0.0006 


+0.0132 
-0.0089 


+25.7 
-18.5 


0.0932 


+0.0077 
-0.0077 


+0.0008 
-0.0008 


+0.016 
-0.011 


+19.6 
-16.0 


950 


0.0457 


+0.009 
-0.007 


+0.0004 
-0.0004 


+0.0084 
-0.0066 


+26.1 
-20.1 


0.0564 


+0.0047 
-0.0047 


+0.0006 
-0.0006 


+0.0104 
-0.0069 


+20.7 
-16.4 


1000 


0.0279 


+0.005 
-0.004 


+0.0003 
-0.0003 


+0.0053 
-0.0036 


+25.9 
-19.3 


0.0348 


+0.0029 
-0.0029 


+0.0004 
-0.0004 


+0.0066 
-0.0043 


+21.2 
-16.6 



(%) 
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In Tables EV- VII we use mw' /2 and 2mw' to estimate the uncertainty due to the choice of scale fj,. While this is 
a reasonable estimate, we show in Fig. H the scale dependence of the LO and NLO cross section for a typical mass 
mw = 500 GeV at run I of the Tevatron over the range fi = mw' /5-5mw'- Even though the scale only enters the 
LO cross section though the PDFs, the NLO cross section is much less scale dependent than the LO cross section. It 
is apparent from this figure that only at unnaturally small scales do the LO and NLO cross sections agree. This is 
typical of single-top-quark production [ ^J39| , and Drell-Yan-like processes in general, where initial-state corrections 
act like new production modes. 




0.2 0.5 1 2 5 

jj/mw' 

FIG. 3. Scale dependence of the leading-order (LO) and next-to-leading order (NLO) cross section pp 
niw' = 500 GeV at run I of the Tevatron {^/S = 1.8 TeV). 



W 



tb; tb for 



In Figs. ^ and ^, we show the LO and NLO single-top-quark cross sections, with W'j^ bosons, at the LHC as a 
function of mw' ■ The cross sections are roughly a factor of 30 larger at the LHC than at the Tevatron due to the 
larger quark luminosity. For W' masses above 1 TeV, the cross section as a function of mass decreases slowly, leading 
to a much larger mass reach at the LHC. For example, the cross sections where quark and lepton decays are allowed 
are 580, 62, 9, and 2 fb, for 2, 3, 4, and 5 TeV W bosons, respectively. Depending on detector performance, limits of 
3-4 TeV on the W' mass should appear within a year or two of running at the initial luminosity of 10 fb~^ per year. 



1000 




mw 



FIG. 4. Leading-order (LO) and next-to-leading order (NLO) cross section pp - 
LHC {^fS = 14 TeV), when the decay to leptons is not allowed. 
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tb, tb as a function of mw' at the 
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B. Parton distribution function uncertainties 



Because we are interested in heavy particles, the parton distribution functions (PDFs) wiU be probed at large values 
of the proton's momentum fraction. Hence, as opposed to Drell-Yan or s-channel single-top-quark production where 
the PDFs are well understood, we expect the PDF uncertainty to play a significant role in the total uncertainty of 
the cross section. In order to estimate the effects of this uncertainty we use a modification [^,0 of the "tolerance 
method" implemented in the CTEQ6 PDFs |||. 

The tolerance method is based on diagonalization of the matrix of second derivatives for ^ (a Hessian matrix) near 
the minimum of for the PDF fits |43|| . Since is approximately parabolic near its minimum Xq, hypersurfaces of 
constant are hyperellipses in the space of the original 20 PDF parameters {a;}. These hyperellipses are transformed 
into hyperspheres by a change of coordinates \ai\ — > {zi}, i = 1,...,20. The tolerance method assumes that all 
acceptable PDF sets correspond to a x^ that does not exceed their minimal value Xo more than by T^. As a result, 
the acceptable PDF sets have {zi} within a sphere of the radius around {^^(xo)} = ^^Vi- I'^ principle, T should 
be chosen to correspond to a Itr deviation of the fit. However, for simplicity we present results for T = 10, as used 
by the CTEQ6M101 - CTEQ6M140 PDF tables given in Ref. ||. 

The PDF uncertainty for an observable O is the maximal change in O as a function of variables {z;} varying within 
the tolerance hypersphere. The CTEQ6 paper estimates the variation of O by using a master formula 



. t ^Ol where SO. . . ^ O^^^t) O^^-^ 



and f = 10 is a step in the space of Zi. Here, 0(z°, . . . , z^' ± . . . , Z20) is denoted as 0(z° ± t). Eq. ( |45| ) is a good 
approximation for the PDF fits, but is less useful for observables, e.g. cross sections. 

The difficulty with Eq. ( ^5| ) is well exemplified for the problem at hand by Fig. |^. Here we see the uncertainty for 
each of the 20 pairs of PDFs as a function of the parameters Zi for mw' = 500 GeV. If we apply Eq. (^) we would 
predict an uncertainty of ±5.9%. However, 0{z^ + ^) ~ ^(z? — t) « for nearly half of the parameters, even though 
Fig. ^ shows that these provide a large deviation. The uncertainty in Eq. (pq ) appears to be more sensitive to Zi 
which have a small absolute effect, but happen to cause a shift in the prediction which changes sign, than to large 
fluctuations. This defect of the master formula is a result of the simple observation that the PDF set that minimizes 
the uncertainty in a given observable O is not necessarily the same as the one that minimizes the fit to the global 
data set. 
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FIG. 6. Uncertainty in the cross section for mw' = 500 GeV for each pair of PDFs using CTEQ6Mlxx. Circles are odd 
PDF sets. Squares are even PDF sets. 



In order to obtain a better estimator of the uncertainty of a generic observable O, we introduce |^,^ the "modified 
tolerance method" (MTM) master formula. We define the maximum positive and negative errors on an observable O 
by 



^(max[O(z0 + t)-O(z«),O(zf-t)-O(z0),0]) , (46) 



g(max[O(z0)-O(z0 + t),O(z0)-O(z«-t),0]) , (47) 
1=1 

where the "tolerance" T is the same scaling parameter that determines the overall range of allowed variation of x^, 
and we use T = < = 10 as in the CTEQ6 PDF sets. In Eqs. (^, we sum the maximum deviations on the observable 
in each of the parameter directions, and hence retain both maximal sensitivity to the parameters that vary the most 
and estimate the range of allowed values of the cross section. Given the case in Fig. ||, we determine that mw' — 500 
GeV has an uncertainty of +8.6—6.2% — half again as large on the high side as estimated by the CTEQ method. 
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C. NLO Differential spectra 

In the course of an experimental analysis it is common practice to normalize the number of events generated in 
an event generator by a "K-factor". This K-factor is generally taken to be the ratio of the NLO to LO inclusive 
cross section. It is then assumed that turning on successive gluon radiation, called showering, in an event generator 
reproduces the correct shapes for the NLO distributions. While this is a reasonable approximation for the soft 
radiation that accompanies the jets, it does not always reproduce the spectrum of additional well-separated hard jets 
|l^ , [45t , or their effect on the primary jets. Since experimental reconstruction efficiencies depend on jet energies, it is 
important to have the most accurate prediction possible. 

There are two main benefits of looking at the fully differential cross section. The first is that this provides a check 
on the jet distributions that come from the event generators. The second is that we see immediately whether the 
kinematic regions where perturbation theory breaks down are relevant to the measurable range of the distributions. 
The shape of the transverse momentum is especially important in studies with &-tags, since the tagging efficiencies 
depend on this variable. 

In Figs. and H we show the transverse momentum (pr) distributions for the top-quark and bottom-jet, respectively, 
for W'^ ^ tb, and mw' = 500 GeV at run I of the Tevatron. For comparison we plot both the NLO distributions 

and the LO distributions times a NLO K-factor of 1.29. The &-jet is reconstructed using a cluster algorithm ]4^ ] 
to provide an infrared-safe way of combining partons. We use a kx cone size of AR = 1, similar to a fixed cone size 
of 0.7. 
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FIG. 7. Transverse momentum prt distribution of the top-quark at NLO, and LO times a K-factor, for mw' = 500 GeV. 
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FIG. 8. Transverse momentum pxb distribution of the 6-jet at NLO, and LO times a K-factor, for mw' = 500 GeV. 



As is generally the case, the pT spectra of the 6-jet and top-quark are somewhat softened at NLO. In the top-quark 
and &-jet pT distributions in Figs. and ^ there is a much improved behavior near the W' resonance at NLO. Aside 
from going to higher orders, it may be possible to further improve the shape in the resonance region by replacing the 
Breit-Wigner with a more dynamical form . While there may be minor deviations at small transverse momentum 
due our choice of a massless bottom quark, the cross section vanishes at low pTb, and 6-jet cannot be measured below 
about 10 GeV, so the effect will not be visible in the experiments. The pseudorapidity distributions of the top-quark 
and 6-jet are shown in Figs. ^ and p^, respectively. The shapes of the distributions are so similar, that the LO curves 
are completely hidden by the NLO curves. 
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-4-3-2-101234 
FIG. 9. Pseudorapidity rjt distribution of the top-quark at NLO, and LO times a K-factor, for mw' = 500 GeV. 




FIG. 10. Pseudorapidity -qt, distribution of the &-jet at NLO, and LO times a K-factor, for mw' = 500 GeV. 

An important consideration is the effect of the choice of cone size on the shape of the transverse momentum 
distributions. In Fig. |ll|, we show the ratio of the pt distributions of the 6-jet for two common choices of Ai? — 0.54 
and 1.35 (which are similar to fixed cones of size 0.4 and 1.0, respectively) to our default choice of 1.0. When comparing 
distributions from event generators with NLO calculations, we must use the same jet definition. Otherwise, there can 
be a systematic shift in the shape of the distributions of 10-20%. This size of this effect is potentially larger than the 
effect of the jet-energy resolution on the reconstruction of jet shapes. 
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FIG. 11. Ratio of distributions of the transverse momentum of the fo-jet using different kr cone sizes at NLO. 



The overall effect of the choice of cone size is only important for W' production right near resonance. The cross 
section is fairly small at lower transverse momentum, and thus changes in the shape in that region are not relevant. 
The effect at low pt is further suppressed in an experiment, because 6-tagging efficiency tends to be smaller at low pT 
as well. In Fig. |l^ we show the effect of different cone sizes on the prt spectrum near the its peak at ~ 215 GeV, which 
corresponds to a 500 GeV W resonance. Experimental analyses will concentrate on this resonant region. Hence, a 
mismatch of jet definitions could cause a systematic mistake in the comparison to the theoretical cross section of as 
much as ±20%. 




While the individual distributions for the W boson are important, in an experimental analysis we want to fit for 
the mass of the W boson. A likely strategy to find the W is to look for a peak in invariant mass distribution Mj^ 
of the top-quark/antibottom-jet pair. We do not perform a full analysis here, but simply show in Fig. ^ the M^i 
distributions for mw' = 500 GeV. The LO cross section is shown times the standard K-factor, but is divided by 2 to 
fit on the figure. While the central value of the peak has not shifted by much, there is a very large tail below the mass 
of the W' boson. Roughly half of the cross section is below the peak predicted at LO. Hard radiation at NLO has a 
much larger effect on the correlated distributions than we might naively expect from the 4% increase in the width. 
Almost 2/3 of the LO distribution falls in a mass window of 500 ± 10 GeV, but only 1/2 of the NLO distribution falls 
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in a mass window twice as wide, 500 ± 20 GeV. Hence, when considering the effectiveness of M^l, the reconstructible 
signal over background may be a factor of 2 smaller than predicted with a LO calculation. 




300 350 400 450 500 550 600 
M^t, (GeV) 

FIG. 13. Invariant mass distribution M^i at NLO, and LO times a K-factor, for mw' ~ 500 GeV. 

The two distinguishing features of the W cross section are that the 6-jet and top-quark are each at a much larger 
transverse momentum than the dominant backgrounds, and there is a peak in the top-quark/antibottom-jet invariant 
mass distribution. While an optimized phenomcnological analysis of the s-channel singie-top-quark production cross 
section at the level of Ref. Q is beyond the scope of this paper, we estimate the reach in W mass using run I 
data in the following way: We begin with the recent limit on the s-channel single-top-quark cross section of ~ 17 pb 
[^,|o|. Assuming that the W boson adds directly to the standard model cross section, this immediately places a 
competitive bound of mw ^ 380-410 GeV (depending on the allowed decay modes). To improve this bound, we note 
that the invariant mass Myyi,^ of the dominant Wbb background (B) scales like l/M^^^^. Hence, by reconstructing the 

W' invariant mass, and scaling the cross section limit by VS, we predict a 95% confidence-level limit of mw' ^ 525- 
550 GeV should be attainable in a dedicated analysis of run I data. In run II of the Tevatron the s-channel single- 
top-quark cross section will be measured to ±20% with 2 fb~^ of integrated luminosity [Q. Given the very low 
backgrounds at large invariant masses, limits of 800-900 GeV may be reachable for W bosons with standard-model 
size couplings. 



IV. CONCLUSIONS 



We present a fully differential calculation of W' boson production and decay at next-to-leading order in QCD. 
We demonstrate that the couplings of the W to fermions factorize through NLO in the width, and in the complete 
two-to-two cross section. Any model that contains a new charged-current gauge particle, with arbitrary vector and 
axial- vector couplings, is described by the analytic results of Sees. || and HI, In particular, the completely differential 
NLO cross section may be used to predict jet distributions for a given model in terms of j. |, Rt, and Ru- In most 
models, left-right mixing is either highly suppressed, or identically zero, which leads to a value of Ru = 0. Hen ce, by 
calculating Rt, we may translate the results for right- or left-handed W bosons, appearing in Sees. HI A and HI C , 
into predictions for these models as well. 

We use our calculation to estimate the effect of a left- or right-handed W boson on the single-top-quark cross 
section. We show that the dominant uncertainties in the theoretical prediction come from the top-quark mass, higher- 
order QCD, and the parton distribution functions. In order to determine the PDF uncertainties, we propose the 
use of an improved modified tolerance method (MTM) that attains maximal sensitivity to the variance of the PDFs, 
and allows us to predict asymmetrical uncertainties. We present the NLO transverse momentum and pseudorapidity 
distributions of the top-quark and bottom-jet. We describe the effect of jet definitions on these distributions, and 
point out the strong sensitivity near the W resonance region. We show that the correlated effect of additional hard 
radiation on the reconstruction of the W' mass peak can be much greater than naively expected from the increase in 
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the W width. 

We predict that the direct search hmit on the mass of the W' boson may be improved by using data from runs I 
and II of the Fermilab Tevatron. By reconstructing the top-quark/bottom-jet invariant mass Mtb in the single-top- 
quark sample, a run I 95% confidcncc-lcvel limit of mw' ^ 525-550 GeV should be attainable for standard-model-like 
couplings. At the end of run II, this limit may be pushed to 800-900 GeV. The signal at the LHC should be large 
enough to reach limits of at least 3-4 TeV. While these estimations are based on the known scaling of the backgrounds, 
a dedicated phenomcnological analysis would be useful. 

We conclude by applying our results to one specific example. In the top-flavor model of Rcf. the left-handed 
coupling in the W'-t-b vertex may be written as 

g L sin C g tan (t>Vtb. (48) 

The coupling to leptons, or the first and second generation of quarks, can be written 

gL sin C Vj'^f^ = 5 cot Vf^^ , (49) 

where g is the standard model SU(2)i coupling, and Vf^f^ is the CKM matrix or diagonal, for quarks or leptons. 



respectively. For masses of around 500 GeV, the angle 4> is restricted to be small |5l|, sin^0 < 0.05. If we choose 
sin^0 — 0.05, the branching ratio into top quarks is 99%; virtually all of the decays of the W will be into top quarks. 
The enhancement in branching fraction is exactly compensated by a decrease in production luminosity. The single- 
top-quark cross section has Rt — tan^cf) cot^^i?^^. Thus, the direct search limit for this top-flavor model will be 
identical to the one derived with standard-model couplings. If we wish to set a limit using a smaller value for sin^c/), 
the W width will be the limiting factor in the mass reconstruction, and we should reevaluate the jet distributions in 
that scenario. By calculating | V^, j, |, _Ri, and and using our next-to-leading order results, we may now accurately 
compare any model containing a W boson to direct experimental searches. 
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